Peculiarities of dynamics of Dirac fermions associated with zero-mass lines 
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Zero-mass lines result in appearance of linear dispersion modes for Dirac fermions. These modes 
play an important role in various physical systems. However, a Dirac fermion may not precisely 
follow a single zero-mass line, due to either tunneling between different lines or centrifugal forces. 
Being shifted from a zero-mass line the Dirac fermion acquires mass which can substantially influence 
its expected "massless" behavior. In the paper we calculate the energy gap caused by the tunneling 
between two zero-mass lines and show that its opening leads to the derealization of linear dispersion 
modes. The adiabatic bending of a zero-mass line gives rise to geometric phases. These are the 
Berry phase, locally associated with a curvature, and a new phase resulting from the mass square 
asymmetry in the vicinity of a zero-mass line. 



I. INTRODUCTION 

Zero energy states are well known for the one- 
dimensional eigenvalue-problem for the Dirac particle 
with a spatially dependent mass. Such a problem natu- 
rally arises in various contexts, both in condensed matter 
and in high-energy physics^. These are super-symmetric 
quantum mechanics^, fractional charge^ and solitons in 
polyacetylene^ 4 . For the monotonous mass distribution 
the wavefunction corresponding to the zero energy state 
is localized in the vicinity of the point were the mass 
vanishes. 

Let us consider the two-dimensional case. Assuming 
that the mass depends on a single variable only, say y, 
we come back to the one-dimensional problem. How- 
ever, in two dimensions the motion is allowed not only 
along the y-axis, but also along the perpendicular to it 
x-axis. If the mass vanishes at the point yo, the line 
(x, j/o ) is the zero-mass line (ZML). Our two-dimensional 
problem reduces to the one-dimensional problem com- 
pletely if we assume that the particle does not move along 
ZML. For non-zero values of the momentum p x along 
ZML the energy of the particle is given by the linear re- 
lation, E = ±p x , i.e. the one-dimensional zero-energy 
state becomes a linear dispersion mode (LDM) in two 
dimensions. 

LDMs naturally appear as edge states for inverted 
band semiconductors^ as well as for a certain model of 
the quantum Hall effeclP. Though single-particle LDMs 
in quantum Hall regime might not be a sufficient de- 
scription for conventional two-dimensional electron gas^, 
this description seems to be adequate for the case of 
narrow graphene ribbons 8 . Another recent examples 
where ZML are essentially involved are given by topolog- 
ically protected edge states in CdTe/HgTe/CdTe topo- 
logical insulators*"^, graphene on boron nitride^^, 
and chemically functionalized graphene^HUl where for 
an inhomogeneous functionalization the mass term can 
in general ch ange i ts sign. LDMs arise in gapped bi- 
layer graphene^HUl and chiral p-wave superconductors^!. 
Changing a width of CdTe/HgTe/CdTe quantum well or 



applying a gate voltage one can also create ZML in bulk. 
Though our consideration will be formally applicable in 
all these cases, for the sake of definiteness we assume 
in this paper that LDMs relate to edge states in two- 
dimensional topological insulators. 

An effe ctive dy namics of charge carriers in topological 
insulator d 10 * 11 * 13 ! is governed by the Dirac Hamiltonian 
with a spatially dependent mass term. This term van- 
ishes along ZML giving rise to LDM. This mode lies in 
a gap for bulk states. Along a single straight ZML the 
linear dispersion mode propagates only in one direction. 
This unidirectional edge mode is very similar to unidirec- 
tional edge states in conventional semiconductors placed 
in the high magnetic field. It is well known, that these 
states support the quantum Hall currentPH The exis- 
tence of unidirectional channels in topological insulators 
in zero magnetic field is referred to, by analogy with th e 
quantum Hall effect, as quantum spin Hall effedP^HSH. 

From the physical point of view a mass square land- 
scape in a topological insulator forms a waveguide around 
ZML (see Fig. [I]). LDM in such a waveguide corresponds 
to zero transversal momentum, i.e. to the rest in the di- 
rection transversal to ZML. This stands in stark contrast 
to the conventional Schrodinger particle in a soft-walls 
waveguide, where the lowest transversal energy is posi- 
tive. 

LDMs in topological insulators are known to be topo- 
logically protected against scattering by non-magnetic 
impurities^. Topological protection results from the 
spatial separation of states traveling in different direc- 
tions, thus the backscattering should be attenuated by 
the probability to tunnel through the bulk of a sample. 
This peculiarity is very similar to the behavior of con- 
ducting states in quantum Hall effect. In ideal system at 
zero temperature it is the tunneling between edge states 
which determines an accuracy of quantum Hall plateaus. 

The assumption about the spatial separation of states 
traveling in different directions does not take into ac- 
count peculiarities of quantum tunneling. Indeed, let us 
consider a ribbon with two parallel ZML at its opposite 
edges and assume that mass does not depend on the vari- 
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FIG. 1: Particles corresponding LDM on the mass square 
surface. Red and green spheres illustrate edge states with 
large longitudinal momenta. They are localized at a single 
zero-mass line. Blue spheres illustrate states with small lon- 
gitudinal momenta. Such states are localized at both ZMLs 
simultaneously. 



able along the line. One can expect that the crossing of 
straight lines corresponding to a linear dispersion turns 
into the avoided crossing. Obviously, this effect does not 
depend on the symmetry: a transition from a crossing to 
an avoided crossing caused by tunneling is generic. 

Due to the gap opening the LDM cone turns into two 
branches, almost linear at large momenta. Let us con- 
sider an upper branch. For large negative momenta the 
wavefunction is localized at one edge of the sample and 
for large positive momenta it is localized at the other 
edge. Since the wavefunction smoothly depends on the 
longitudinal momentum, we should conclude that at a 
certain longitudinal momentum amplitudes of the wave- 
function are comparable at both ZMLs. At this point 
states traveling in different directions are not spatially 
separated and the topological protection may break (see 
Fig. [I]). In Section III we illustrate this effect by an ex- 



ample of a symmetric mass distribution. 

The tunneling between ZMLs has been studied 
experimentally 2 -. In this work the differential con- 
ductance measurements in the integer quantum Hall 
regime between two parallel edge channels were reported. 
Theoretically^ one considered the tunneling between 
Quantum Hall edge states via the Landau-Zener-like 
mechanism: it was assumed that the main contribution 
comes from the vicinity of a point were edge channels are 
the most close to each other. 

Topological protection of a single ZML should lead to 
a substantial change of dynamical properties of LDM. 
Indeed, let us consider an adiabatically bent ZML. For 
a conventional Schrodinger particle the bending implies 
a geometric potential, propo rtional to the square of the 
curvature at every pomt r 8 l 29 l This barrier results in the 
appearance of a trapped mode and backscattering. Both 
effects are forbidden for LDM due to the topological pro- 
tection. In the Section |IV| we show that the bending of 
ZML, indeed, tends only to the emergence of geometric 
phases. 



II. LINEAR DISPERSION MODES AND 
ZERO-MASS LINES 

We assume that the dynamics of charge carriers is gov- 
erned by the Hamiltonian 



H = <T x p x + (T y py + a z m(y), 



(1) 



where the mass m(y) vanishes along the line y = 0. Here 
we put H = v = 1, where v is the Fermi velocity. For a 
given energy E the wavefunction 'F obeys the equation 



= £<F. We write * 



x(y), 



WxPx + <J y Py + <J z m{y)]x{y) = E\(y). (2) 
In a matrix form |2| reads 

m - E p x - d y ) ( \ 



Px + d y -m - E 



X2 



= 0. 



(3) 



Let us now sum up equations in ^ and subtract the first 
one from the second. Then we find 



p x - E dy + m 
d y -m Px + E 



0. 



(4) 



where we introduced the notations rji = (xi + X2)/v2, 
V2 = (xi —X2)/V% (see Appendix |a| for details). We can 
reduce Q to scalar Schrodinger equations^ 



-dy + m(y) 2 + m'(y)]r)i = A771, 



(5) 



[-d 2 y + m{y) 2 - m' (y)] m = Ary 2 , (6) 

where E 2 = p 2 + A. Functions 771 and 772 are not inde- 
pendent. Connection formulas read 

(E + p x )r) 2 = (rn- d y )rji, (7) 
(E - p x )i]i = (to + d y )rj2. 

Equations |5])-(|6| can be written in the form 

(to + d y )(m - d y )r)i = Xr)i, (8) 
(to - d y )(m + dy)r) 2 = Xi] 2 . 

Multiplying the first equation in pi by ryi , the second by 
1J2 and integrating over y we find 



a|MI 2 = |I(™- WI 2 >o, 
\\\m\\ 2 = \\(m + d y ) m \\ 2 >o, 



(9) 



whence A > 0, since 771 and r\2 can not vanish simultane- 
ously. In ([9]) we denoted 



vf(y)dy. 



(10) 



Let us consider a case when m monotonously depends 
on y, say m'(y) > 0. Then equations Q comprise a 
LDM, for which E = —p x . It is clear that 771 should be 
identically zero, since ^ can not have a zero eigenvalue 
if the potential is always positive. From Q for 7/1 = 
we have E — —p x and 



772(2/) = exp ( - J dy'm(y) 



(11) 



FIG. 2: Tilted double well potentials. Solid and dashed lines 
correspond to Vi, V2 respectively. An existence of the up- 
turned well can lead to the Arago effect: focusing in a shadow 
region 



FIG. 3: The schematic view of the dispersion relation and 
the corresponding wavefunction reconstruction (insets). One 
sees that at zero longitudinal momenta the wavefunction is 
not localized at a single ZML. 



III. TWO ZERO-MASS LINES 



Let us consider a mass distribution 



™>{y) =y 2 - a 2 . 



(12) 



It mimics a sample with two parallel ZMLs at y — ±a. 
The infinitely large values m(y) in the limit y ±oo 
can be treated as two edges of the sample. Effective 
potentials 



vi(y) = m(y) 2 + m'(y) = (y 2 - a 2 ) 2 + 2y, 

v 2 (y) = m(y) 2 - m '(y) = (y 2 - fl2 ) 2 - 2 y, 



(13) 
(14) 



entering equations (|5|, in this case correspond to 
tilted double wells (see Fig. [2). 

It is easy to see that A = does not belong to the spec- 
trum of Eq. ([6]). Indeed, the exact solution (11 1 exponen- 
tially decays when y — > oo and exponentially grows when 
y — > — oo. Due to the conservation of Wronskian another 
linear independent solution of ^ exponentially grows 
when y — > oo and exponentially decays when y — > — oo. 
Similarly we prove that A = does not belong to the 
spectrum of Eq. 

Nonzero A leads to avoided crossing of branches E = 
±p x , corresponding to well separated ZMLs (se e Fig. [3]) . 
The dispersion relation reads E^(p x ) = ±yjp x + A n . 
The splitting between branches E^{p x ) is equal to 2\/Ao, 
where Ao is the lowest eigenvalue of the Schrodinger equa- 
tions ([5| , ^ . We put p x = and take the branch corre- 
sponding to Eq (0) = p y = V^o ■ Equations ^ give 



d v )m(y) 

dy)Vi(y) 



(15) 



where fji(y), fj 2 (y) correspond to p x = 0. If fj^y), fj 2 (y) 
are known then the functions 771 (y) , 772 (y) for p x ^ can 
be reconstructed as 

vx(v) = ( E + Px)fji(y), m(y) =p y m(y)- (16) 



Let us multiply the first equation in ( |15[ ) by f}i(y), the 
second one by fj 2 (y) and subtract the second result from 
the first one. We obtain d y [f]i{y)f] 2 {y)] = P y [fj 2 (y) — 
fj 2 (y)]- Integrating the last equality over y from minus to 
plus infinity we find ||^i|| 2 = ||^2|| 2 , which suggests that 
Xi,2(y) = [ ? 7i(y)± r yi(y)]/v / 2 are not localized at a certain 
ZML at zero longitudinal momentum p x , but rather have 
comparable amplitudes at both ZMLs. In contrast to the 
double well problem (see® , p. 183) this effect is not a 
consequence of a spatial symmetry. The derealization at 
zero longitudinal momentum may destroy the topological 
protection against disorder. 

Let us now look how the localization appears at non- 
zero longitudinal momenta. For \p x \ 3> p y we find E + 
Px Py if Px > and E + p x <C p y if p x < 0. Thus for 
large positive longitudinal momenta X1.2 are localized at 
one ZML and for large negative longitudinal momenta 
these functions are localized at another ZML. 



For symmetric mass distribution (12) we can obtain 



an analytic estimation for the spitting. Let us change y 
to — y in (15). Taking into account that m{y) is an even 



function we find 

(m - d y )fj 2 (-y) =p v m(-y), 
{m + dy)fji(-y) =p v m(-y). 

Thus we conclude fji(—y) = fj 2 {y)- Let us multiply the 
first equation in (15) by fji(y) and integrate over y. We 
obtain 



Py\\m\ 



dy f] 2 {-y){m + d y )fj 2 {y) 



— OO 
OO 



= 2 / dyf, 2 (-y){m + d y )fi 2 (y)+fi 2 2 (0) 
10 



2 Py dyvf(y)+ij 2 (0) 



(17) 



Since the function 771 obeys the Schrodinger equation ^ 
with potential ( 13 ) for small p y it remains exponentially 
small for any y > except the vicinity of the point y = a. 
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In this vicinity fji grows due to the focusing in a shadow 
region (an effect similar to the Arago spot, Fig{2|, but 
nevertheless remains small. Therefore we can neglect an 
integral in the last equality in (17) and write 



where ( 1 1 



Py\\m\\ 2 ^v 2 2 (0), (18) 

was used to approximate 772 in the region y > 



0. Using (11) again we hnd 



Py 



r°o / rv 

J dy cxp ( -2 J dy'm{y') 



= 1. 



(19) 



method. It gives 
<7r 



The integral in (|19| can be computed using the Laplace 

dy'miy')) = 1, (20) 



whence 

Py = y^ cx P ( 2 J dy'm{y')*j = y^cxp (-4a 3 /3) . 

(21) 

Thus for the splitting we obtain 



2\p y \ = 2^/- exp I 



-4a 3 /3) 



(22) 



Though the explicit estimation ([22]) of the splitting 
holds only for the specific symmetric mass distribution 
(12), one can expect that for a generic case the splitting 



in LDM for a system with two ZMLs is proportional to 



exp 



\m(y)\dy , 



(23) 



where a\ < ai are positions of ZMLs. In contrast to the 
famous double-well problem, the ground state A > of 
([5]) , ([6]) is still determined by the tunneling, since for any 
m(y), every well treated apart generates a zero eigenen- 
ergy. For the double well the semiclassical degeneracy 
persists in the symmetric situation only. 



IV. BENT ZERO-MASS LINE 

Let us consider a bent ZML given by {x, y} = R(r), 
where t is a natural parameter, i.e. |R'(t)| = 1. In the 
vicinity of this line we introduce new variables r, £ by 
the equality {x,y} = R(r) + £n(r), where n is a unit 
normal vector on the curve at the point r. In curvilinear 
coordinates |l]) reads 



H 



io-R'(r) d 

l-ik{T)dr 



j(7n(r) 



d_ 



a z m, 



(24) 



where k(r) = — (R',n') is the curvature at the point r 
and arh = a x b x + <J y b y for a vector b = {b x ,b y }. Since 
the Jacobian 



J 



D(x,y) 



1 - fc(r)£ 



(25) 



is not unity, we introduce a new wavefunction 

* = v/l-fc(r)£*, (26) 

which in curvilinear coordinates has a "conventional" 
normalization condition: 



/ dTd£{¥W) = 1. 

Jv 



(27) 



Then the stationary Dirac equation reads = E^> , 

~ (TR'(r) „ , , d 

H = - ^wPr ~ «crn(r)— + a z m 

ik(rn(T) icrR'(r)£fc'(r) 



2(l-£fc(r)) 2(l-£fc(r))2- 
In the case \m' T \ <C |m4| we find H ~ Hq + Hi, 

d 

H = ffR'(r)p T - icrn(T)— + a z m, 

Hi=<T*!{T)£,k{T)p T -~<Tn{T). 

In the adiabatic approximation the effective dynamics 
is one-dimensional along ZML. It is governed by the ef- 
fective scalar Hamiltonian L ~ Lq + Li. The symboP^ 
Lo(pt,t) of Lo is an eigenvalue of the problem^! 



(28) 
(29) 



d \ 
(tR'(t)p t -i<rn(r)— +a z m \ x(Pt,t) 

= L {p t ,t)x(Pt,t). 



(30) 



Using the notations n = {711,712}, R' = {712,— «i} we 
obtain 



m — La p T — c\ 
p T + —m — L 



Xi 
Xi 



0, 



(31) 



where xi — ( n 2 — ifii)xi- After to the replacement r 
x, £ — > y (31 ) coincides with (Q. From the expressio" 3 ' 



{X &p T drj 



we find 

Li = 



2%37 + ^»>e*(r)ft- " 



(32) 



(33) 



Here (-) g m eans the integration over £. To pass from 
(32) to (33) we have chosen x = {XI1X2} to be a real 
function and used the following equalities: x^( <Tn )x = 
0, (xh'pji = 0, xH<r&)x = 2 X ix 2 , X H^)drX = 
(d T -ik)xiX2, (XiXa)s = Pr/2L , (£o)p T =Pr/L a . 

The solution ii of the effective longitudinal equation 
L-0 = £:t/i read^S 



Lo 



(34) 



= / p T dr - 



^o(xi£x 2 )« - ^ fc(r)dr. (35) 
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FIG. 4: Illustration of the closed ZML. Solid arrows show 
the tangent vectors, dashed arrows correspond to x-axis. 



It relates to the solution of the Dirac equation as ^> — 
Xip- From the definition of curvature we have fc(r)rfr = 
— (R',dn). Introducing the angle <j> between R' and x- 
axis we find from the last equality fc(r)dr = dcj). This 
gives 



= / p T dr- J Lo(xiCX2)efc(r)dr 



Acj) 
~2~' 



(36) 



where Acf> is the total rotation of the tangent vector. 

From ( 34 ) it is easy to find the quantization rule. For a 
closed ZML without turning points (see Fig. [4| we obtain 



1 

2^ 



p T dr 



E 
2tt" 



w 



(37) 



where p T = ±\/ E 2 — A(t), A(t) = L 2 (p T — 0, r) and w 
is the winding number indicating how many times the 
tangent vector turns around a fixed point. For a curve 
without intersections w = 1. 

For LDM L = E = —p T provided that > and 
(33) gives 



L i = (XitX2)tk(T)p T 



(38) 



The first term in ( 38 1 describes the energy gain due to 



the displacement of the Dirac fermion from ZML caused 
by the confinement asymmetry. This energy gain is the 
consequence of the centrifugal force. One sees that in 
the considered approximation the centrifugal force itself 
does not shift the particle since its energy is assumed to 
be relatively small. For a higher energy one can use the 
approach developed irPSl 



The term k/2 in (38) is a geometric poten tial associ- 
ated with a curvature. It is well-knowr l 28 * 29 ! that for a 
Schrodinger particle in a bent waveguide the geometric 
potential is negative and proportional to k 2 . It always 
leads to the formation of curvature induced bound states. 
On the contrary, for LDM the geometric potential results 
in the appearance of a geometric phase only. Indeed, the 
effective longitudinal wavefunction ip p4| reads 



= -Et - E 



rl>(r) = 
(Xi£X2)?fc(T)dr + 



Ac/> 
~2~' 



(39) 



The absence of trapped states for the Dirac fer mion c an 
be seen as a manifistation of the Klein tunnelin g ^ 34 * 35 ! for 
the massless LDM. A similar effect was already found 
irP^for a mode with a liner dispersion in a bent graphene 
ribbon. 



For the LDM quantization condition (37) can be sim- 
plified to give 

El E f , _ _ . , . , , w . 

-^-^ f(x^X2) i k(r)dT = n--. (40) 



Here I is the length of the closed ZML. From (40) we 
obtain the semiclassical spectral series E n : 



E n — — - 



(MX2> ? fc(r)dT. (41) 



Half-integer numbers in the first term of (41 ) can be seen 



as a manifestation of the Berry phase for the massless 
Dirac fermion. Indeed, along ZML a particle described 
by LDM remains massless, therefore the conventional ar- 
gument that such a particle acquires the Berry phase 
equal to irw can be repeated. 

Using the results of Section |TT] for a single ZML we 
find Xi — ~X2 — %/\/2) where 772 differs from (11) by a 
normalization factor 



N(t) 



This gives 



exp -2 / d?m(T,g) d£ 



(Mxi)i = - J £ cxp f -2 JT d£'m(T, J «■ 

(43) 

More precisely, the integration in ( 42 ) , ( 43 1 should be 
performed between finite limits lying in a sufficiently 
large vicinity of ZML which, on the other side, does not 
contain any other ZML. We completely neglected tunnel- 
ing effects, thus the constructed LDM should be consid- 
ered as an asymptotic of the eigenfunction of the Dirac 
equation in the given vicinity of ZML. The non-zero value 
of (43) means that the Dirac particle is slightly shifted 



-1/2 



(42) 



from ZML due to the local asymmetry of the confine- 
ment. Indeed, if m is an odd function of £ the expression 



(43) vanishes. This happens since the confining potential 
m 4 {T, £) in this case is symmetric with respect to ZML. 

V. CONCLUSION 

In the paper we studied the LDM dynamics of charge 
carriers in topological insulators. We have shown that 
the wavefunction of a charge carrier is localized along a 
single ZML only at large longitudinal momenta. At small 
longitudinal momenta the wavefunction has comparable 
amplitudes at both edges of the sample, which may affect 
the topological protection. We found that the curvature 
of a bent ZML forms a geometric potential, which how- 
ever does not lead to an appearance of trapped modes 
due to the Klein tunneling. 
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7T + i(J x [TTy , 7T 2 J + lOy [TT Z , TT X \ + 1(7 z [K x , TT y 



(Al) 



valid for any non-commuting operators w x , Tf y and tt z . 
We find 

\pl +pl+ m{yf + a x m' (y)} X = E 2 X - (A2) 



Appendix A: On the effective Schrodinger equation 

In this abstract we would like to clarify the origin of 1 
the transformation performed to pass from ^ to Q. ./2 
Let us write equation (pi) in the form <TTtx(y) = Ex(y), 

where tt — {p x ,p y ,m(y)\. We can square this equation and obtain ([5]), ([6]). 



The last equation comprises the only matrix a x , which 
does not depend on y. We write 











" v/2 ( 


-0 



V2 



(A3) 
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